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Abstract: 

The statistics of g-oscillators, quons and to some extent, of anyons are studied and the 
basic differences among these objects are pointed out. In particular, the statistical distri- 
butions for different bosonic and fermionic g-oscillators are found for their corresponding 
Fock space representations in the case when the hamiltonian is identified with the num- 
ber operator. In this case and for nonrelativistic particles, the single-particle temperature 
Green function is defined with g-deformed periodicity conditions. The equations of state 
for nonrelativistic and ultrarelativistic bosonic g-gases in an arbitrary space dimension 
are found near Bose statistics, as well as the one for an anyonic gas near Bose and Fermi 
statistics. The first corrections to the second virial coefficients are also evaluated. The 
phenomenon of Bose- Einstein condensation in the g-deformed gases is also discussed. 
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1 Introduction. 



In the last few years the interest in particles obeying statistics different from Bose or 
Fermi has increased considerably. The observation in 1953 by Green that different 
kinds of statistics are allowed within the context of quantum field theory (QFT) has 
stimulated an extensive study of the so-called paraBose and paraFermi statistical fields |0] . 
The commutation relations for these fields are trilinear in the creation and annihilation 
operators and are characterized by an integer p, the order of the parastatistics, which 
corresponds to the number of particles in a given symmetric or antisymmetric state. 
(If p = 1, one recovers the usual Bose and Fermi statistics). The case when p is not 
an integer has been recently also studied in order to provide theories in which the 
Pauli exclusion principle and/or Bose statistics can be slightly violated. However, the 
corresponding QFTs, which describe such particles, turn out to have negative norm states 
and consequently are not physically acceptable. More recently, another possibility has 
been explored, which corresponds to the case where no assumption is made about the 
parameter p [0-0. The particles which obey this type of statistics are called "quons" 
and the statistics is referred to as " infinite statistics" . They can be realized by different 
commutation relations like the ones of g-deformed algebras. 

On the other hand quantum groups [0 are a subject of great activity at present and 
although their direct physical interpretation is still lacking, it is of particular importance 
to study the possible physical implications of these deformations. These structures which 
first emerged in connection with the quantum inverse scattering theory and solvable 
statistical mechanical models , have found recent important applications in many prob- 
lems of physical and mathematical interest, such as rational conformal field theories |T^, 
non- commutative geometry [0, knot theory |jl2[, quantum superalgebras [0 and so on. 
The concept of g-deformed oscillators (the so-called g-oscillators), derived from the 



contraction of quantum algebra [|T3|], has also been introduced. These g-oscillators can 
formally be defined in any space dimension but then they violate the fundamental axioms 
of QFT on the relation between spin and statistics. In two dimensions, however, a new 
kind of statistics interpolating between bosons and fermions can (without violating any 
of the axioms) exist, due to the appearance of the braid group instead of the permutation 
group when identical particles are exchanged. The particles which obey such so-called 
fractional statistics are called anyons |T^. Such particles have attracted considerable in- 
terest since they appear in realistic systems and in particular, are connected with the 
interpretation of the fractional quantum Hall effect [|TB| and may have relevance in other 



quasi-two-dimensional systems such as high Tc superconductors |jT^ 



Since anyons are two-dimensional particles acquiring a phase upon circling each other 
and returning to their original configuration, their wavefunctions as well as the operators 
that create or annihilate them must be multi-valuated functions of position. This fact 
makes the construction of such operators rather involved. This problem was originally 



studied for field theories in the continuum [18|,|19] but there the precise definition of 



multi-valued operators encounters severe difficulties. The same problem, however, has 



been lately examined for field theories on a lattice where the multi- valued 

anyon operators can be indeed rigorously defined. 

In this paper we would like to study the interrelation between the statistics of q- 
oscillators, quons and, partially, of anyons. At the same time we point out the basic 
differences among the statistics of the above mentioned objects. Along this line we hope 
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to be able to use the techniques, such as the Green function method developed here, 
the Fock space representation, etc., for the g-oscillators, further in the case of fractional 
statistics. 

In sect. 2, the g-deformed bosonic and fermionic oscillators are defined and their Fock 
space representations are given for different forms of their commutation relations. We 
emphasize the basic differences between g-oscillators, quons and anyons in general. Sect. 
3 is devoted to the study of the statistics of g-oscillators. In particular, the statistical 
"distributions" for g-bosons and g-fermions are found in the case when the hamiltonian 
of the system is identified with the number operator. Comments on the results previ- 
ously obtained in the literature, are presented. In sect. 4 the g-deformed single-particle 
temperature Green function is defined and its "periodicity" properties are studied in or- 
der to get an appropriate Fourier decomposition in the Euclidean "time" variable and 
find the corresponding frequency values. The method for evaluating frequency sums is 
also discussed and in particular, is used to calculate the statistical "distribution" for q- 
particles. In sect. 5 we consider an ideal g-gas described by a general noninteracting 
hamiltonian and find perturbatively the equation of state near Bose statistics in different 
cases, namely, nonrelativistic and ultrarelativistic g-boson gas in an arbitrary space di- 
mension. The low-density regime is also considered to find the virial expansion and the 
corrections to the second virial coefficient. Sect. 6 is devoted to the study of a gas of 
anyons near Bose and Fermi statistics. We find perturbatively, in the first order of the 
statistics determining parameter, the equation of state of the anyonic gas by means of 
the temperature Green function approach. Finally, concluding remarks and directions for 
further studies are given in sect. 7. 



2 g-oscillator algebras. 

In this section we shall review briefly the definition of g-deformed bosonic and fermionic 
oscillators and their Fock space representations. As in the classical case, the g-oscillators 



can be obtained from the SU{2)q algebra by a contraction procedure [|T3 



We start by defining the g-bosonic oscillator algebra through the commutation rela- 
tions p|,p^ 



+ + ~N 

aa — qa a = q 



[N,a] = -a, [N,a+] = a+ , (2.1) 

where g G C ; a, a+ and are the annihilation, creation and number operators respec- 
tively. One can construct the representation of (2.1) in the Fock space spanned by the 
orthonormalized eigenstates |n > of the operator A^, 



\n >-- 



|0 > , a\0 >= , A^|n >= n\n > 



(2.2) 



where [n]\ = [n][n — 1]...[1], [0]! = 1 and [n] denotes 



n 



g" - g-" 



g - g 



-1 



(2.3) 
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In this Fock space it is easy to prove that the following relations hold: 

a+a = [N], aa+ = [A^ + 1] . (2.4) 

It is interesting to note that if g = e^*'^/™, then [m] = (see definition (2.3)). Conse- 
quently, the Fock space breaks up into m-dimensional subspaces not connected by the 
operators a and a~^. Each subspace carries an m-dimensional representation of the al- 
gebra (2.1) and can be considered separately. (An example is given below, eq. (3.6)). 
A corresponding phenomenon happens whenever q = exp(2r7r), where r is rational, and 
different from an integer. 

These exist different forms of relations (2.1). In particular, when q is real, it is possible 



to define the operators |22 



JV ,1 , N 



A = q~a , = a+g~ , (2.5) 
which satisfy the commutation relation 

AA+ - q^A+A = 1 . (2.6) 

In terms of the new operators (2.5), the Fock representation (2.2) reads as 

\n >= y ' 10 > , AIO >= , N\n >= n\n > , (2.7) 

where 

[nf ^ q-~'[n] = . (2.8) 

(B standing for bosons). 

Finally in this Fock space we have 

A+A = [Nf , AA+ =[N+ if . (2.9) 

One can also express (2.1) in terms of the usual undeformed Bose operators b, 6^ by means 



of the change |22 



Then we obtain [b,b~^] = 1, [N,b] = —b, [N,b'+] = b+, N = b'+b, i.e., the usual bosonic 
oscillator algebra. 

Let us now introduce the fermionic g-oscillator algebra. The annihilation and creation 
operators of the fermionic g-oscillator /, /"*" are postulated to satisfy the commutation 
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+ = , (2.11) 

[N, f] = -f, [N, /+] = /+ , 

where is the g-fermion number operator. The orthonormalized eigenstates of are 
defined by 

\n >= ^^p^|0 > , /|0 >= , N\n >= n\n > , (2.12) 



n 



where 



and we have 

r/ = [iv]^ fr = [N+iY . (2.14) 

For generic g, this representation is infinite-dimensional. Like in the g-bosonic case, for 
special values of q the Fock space breaks up into disjoint subspaces each carrying a finite- 
dimensional representation of (2.11). If, in particular q = e.^™!^ ^ we have \lvn\^ = for m 
odd, whereas for m = 4k, [mY = and for m = 4fc + 2, k > 1 [^Y = 0; the dimensions 
of the disjoint subspaces being the smallest integer n for which [nY = 0. (The cases 
q = e^""^/^ are special, with e.g. [nY — — i^^'^'^^'^n; they thus give infinite-dimensional 
representations). For q = 1 the Fock space breaks up into two-dimensional subspaces, 
and the Pauli exclusion principle follows from = (Z"*")^ = 0. 
Let us introduce now the modified operators 

F = / , F+ = /+gT , (2.15) 

for q real, which satisfy the commutation relation 

FF+ + q'^F+F = 1 . (2.16) 

Using (2.12) one can build the Fock space representation in terms of these new operators. 
We obtain 

\n >= V ^ |0 > , F|0 >= , N\n >= n\n > , (2.17) 



with 

[nf ^ q--'[nY = ^ \ ^ . (2.18) 



1 - (-1) V" 

g2 + l 



It is worthwhile to remark that in the case of fermionic g-oscillators defined through the 
commutation relations (2.11), there does not exist a change of operators, analogous to 
(2.10), which would allow us to express these oscillators in terms of the usual (undeformed) 
Fermi ones. One can use instead another definition for fermionic g-oscillators |T3|, 



cc^ + qc^c = q^ 



[M,c] = -c , [M,c+] = c+ . (2.19) 
In this case the change of operators 



M 



C = q--c , C+ = c+g~- (2.20) 
would lead to the usual Fermi commutation relations 

CC+ + C+C = 1 , 

and consequently, the only non- vanishing eigenstates of the number operator M are |0 > 
and |1 >= C~^\0 >. (Notice that this conclusion derived from (2.19) is equivalent to the 
relations = = and thus, there is no need to assume them). This means that 
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fermionic g-oscillators defined by the algebra (2.19) have only the trivial (undeformed) 
Fock space representation and thus, they behave as usual fermions. 

To conclude this section, some remarks about systems with several degrees of freedom 
are in order. In the case of bosonic and fermionic g-oscillator algebras (2.1), (2.5) and 
(2.11), (2.15) respectively, the multimode generalization is straightforward. For g-bosons 
we have the commutation relations 

^iaj - ((? - l)'^ij + l)a^«i = Sijq~^' , 
[Ni,aj] = -Sijaj , [Ni,aj] = 6ijaj ; (2.21) 



AA+ - ((g2 _ i)S^^ + i)AjAi = % , 

[N„ A,] = -6,,Aj , [N„ A+] = 6,,A+ ; (2.22) 
while for g-fermions they read as 

/,/; + ((g-l)5., + l)/+/, = <5,,g-^» , 

[N^, /,] = -5i,fj , [iVi, /+] = 5.,/; ; (2.23) 



F,F; + ((g2 - i)S^^ + 1)F+F, = 5,, , 

[iV„ F,] = , [iV„ F;] = 6,,F; . (2.24) 

Recently, an alternative deformed commutation relation, called quon algebra, has been 
considered It is defined by 

aiaj^ — qa^ai = 6ij , (2.25) 

and interpolates between Bose and Fermi algebras as q goes from 1 to -1 on the real axis. 
However, let us note that there exists a distinct difference between the g-oscillators defined 
in (2.22), (2.24) and the quons. In fact, in the case of bosonic (fermionic) g-oscillators 
different modes [i ^ j) commute (anticommute), while quons "g-mute" i.e. satisfy the 
relation given in (2.25). Moreover, in the case of quon algebra no commutation rule can be 
imposed on aa and a^a'^. In effect, relations such as a^aj — qajai = or afoj—qa^af = 
would hold only when = 1, i.e., only in the Bose and Fermi cases. 

As far as anyons are concerned, their braiding properties make the construction of an 
oscillator algebra a difficult and still unsolved problem. In a recent paper the construction 
of anyonic oscillators on a lattice as well as their relation with quantum groups have been 



studied I^T]]. We will not discuss them further here, but it is worthwhile to stress that g- 
oscillators and anyonic oscillators differ substantially from each other. First of all, anyons 
arise only in two dimensions, while g-oscillators can be defined in any space dimension. 
Secondly, as we shall see in section 4, the g-oscillators can be interpreted as the Fourier 
components of local field operators, whereas anyons are non-local objects as a consequence 
of their braiding properties and the essential differences in their commutation relations. 
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3 (/-deformed statistics. 



In order to study the properties of g-oscillators, let us consider now the statistical averages 
and calculate the deformed "distributions" which follow from the g-algebras defined in 
the previous section. 

As is well known the thermodynamic properties are determined by the partition func- 
tion Z, which in the canonical ensemble is defined by 

Z = Tr(e-^^), (3.1) 

where f3 = l/{kT), H is the hamiltonian and the trace must be taken over a complete set 
of states. For any operator O, the ensemble average is then obtained with the prescription 

<0>= ^Trie-f^'O) 

and it is crucial that the cyclicity of the trace be consistent with the algebraic structure. 
In the remainder of this section we restrict ourselves to systems described by a single 
oscillator mode. Multimode systems are treated in sections 4 and 5. 

Two types of hamiltonians are usually dealt with in the literature. The first one is 
identified with the number operator and reads as [Q, 

H = ujN, (3.2) 

(hereafter h = 1), while the second one is defined in terms of the creation and annihilation 



operators, following the classical realization of the harmonic oscillator 

H = ^{a^a + aa+) . (3.3) 
For the g-fermion oscillator the latter can be taken as 

^ = |(/v-/n ■ (3.4) 

In section 5 we shall generalize (3.3) to a general linear combination of a~^a and aa"*". 

The partition function, and hence the thermodynamics of the system, is determined 
uniquely by the spectrum of the hamiltonian. With the choice (3.2), the partition function 
will be that of the harmonic oscillator. 



E e-'- = ^ (3.5) 



n=0 



irrespective of what deformation, Bose or Fermi, we want to choose. Of course, for those 
special values of q for which the Fock space separates into disjoint subspaces each carrying 
a finite-dimensional representation, it is perhaps physically more meaningful to consider 
the truncation to one of these subspaces only; i.e. picking only those terms in the sum 
(3.5) corresponding to the subspace chosen. As an example, when q = e*'^/^, the Fock 
space states |3 >, |4 >, |5 > span a representation R of (2.1) in which 



aR = \ I \ , a^ = \ I \ , Nr= 4 (3.6) 
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(note that a"*" is not the conjugate of a in this representation). The corresponding trun- 
cated partition function is 



^-W^ + e"^^^ + e"^'^" . (3.7) 



Summing over all subspaces reproduces the full partition function (3.5). 

While the thermodynamics for a system with hamiltonian (3.2) is independent of the 
deformation, Green functions like < a^a > will depend on the deformation. For the 
g-bosonic algebra (2.1) we obtain 



< a+a >=< [N] >= ^ENe"^"" 



(3.^ 



< A^A >=< [Nf >= -^^^ . (3.9) 



^2Pu _ (g + g-l)e^'^ + 1 ' 

while for the operators (2.5) obeying (2.6) we have the simple expression 

1 

Similarly, for the g-fermion algebras (2.11) and (2.16) we obtain 



<tf> = < [NY > = (3.10) 

and 

< F+F > = < [Nf > = -T^—, . (3.11) 
The average of is of course always given by the usual Bose-Einstein formula 

n=<N>=^ (3.12) 

(when the average is taken over the complete Fock space), to which (3.8) and (3.9) reduce 
in the case g = 1, and (3.11) in the case g = ± as they should. When g = 1 in (3.10) 
or (3.11), the Fermi-Dirac distribution is recovered. 

When g 7^ 1, the temperature Green functions < a^a >, etc. do not have a direct 
relation to the thermodynamic quantities of the system. In a multimode context, such as 
will be considered in section 4 and 5, this can be understood by noting that the average 
< a\ak> does not describe the average occupancy of the kth level, < N^ >. 

When the hamiltonian is chosen as in (3.3) or (3.4), the partition function cannot any 
longer be computed in closed form for generic g. This case has been investigated in refs. 
121-121. 



4 Temperature Green function approach. 

In this section we shall define the single-particle temperature Green function for g-particles 
satisfying the algebras (2.6) or (2.16). This would mean that the fundamental field oper- 
ators of the theory are the ones connected to the operators appearing in (2.6) and (2.16) 
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(see (4.4)). In the case of the algebras (2.1) and (2.11) the construction of this function 
is a more comphcated task due to the nonhnearity of the commutation relations with 
respect to the annihilation and creation operators, which makes it less obvious how to 
define the Green function. 

Let us start by defining the single-particle Green function as 

Ga/3ix, r; x', r') = - < T,M^, T)V^(f' , r') > , (4.1) 

where the ordering operator Tq orders the operators according to their values of r, 

+ /XX i il'(x,T)'ih~^(x',T') if r > r' 

T,M., rmx'. r')) = I ,;;.'(,-;;;,)V(f , if . < ^ (*-^) 

In this section in order to simplify the notation and consider both g'-bosonic algebra (2.6) 
and g-fermionic algebra (2.16) simultaneously, we shall use q instead of ±g^. Note that 
when q = 1 {q = —1) eqs. (4.1) and (4.2) are the usual definitions for bosons (fermions). 
The indices a and (3 label the components of the field operators and in what follows they 
will be omitted. Finally, the operators 4'{x, r) and ip~^{x, r) are the Heisenberg operators, 
related to the Schrodinger ones %l){x) and ip^{x) through the relations 

V'(f,T) = e^^i/j{x)e-^^ , 

^+{x, t) = e^^7jj+{x)e-^^ . (4.3) 

As a simple example, we shall compute the Green function for a spin-zero noninteracting 
system. In this case the states are plane waves so that 

m^^i:e^''A,, (4.4) 

and the operators A^^A^ are supposed to satisfy the g-oscillator algebra 

AkA+ - qA+Ak = 1 , (4.5) 

[Ni,Ak] = -SikAk, [Ni,At] = 5ikA+ . 
The free hamiltonian will be taken as 

H = J2ekNk, (4.6) 

k 

SO that the following commutation relations hold 

[H, A,] = - Gfe Afe , [H, At] =efe At . (4.7) 
The Heisenberg field operators (4.3) are then 

^(S,T)^^Y.^'''~^''Ak, (4.8) 
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k 

Now, from the definition of the Green function (4.1) and using (4.8) we have 
G°{x,r; f ,t) = -l-J2e'^^^-^'^'^e-^'^+^^''' 

kk' 

x{9{t - t') < AkA+ > +e{T' - T)q < A+Ak >} 
= -1 ^e^*^(*-^')e-^^("-"'){^(r - r') < A^A^ > 

k 

+e{T' -T)q<AtAu>} . (4.9) 

The last equahty in (4.9) follows from the orthogonality of states with different occupation 
numbers {uk}- For the ensemble averages we have according to eq. (3.9) 



< AtAu >= , (4.10) 

and from eq. (4.5), 

< A^At >^l + q< A+Ak >= . (4.11) 



Let us now turn our attention to the "periodicity" properties of the Green function (4.1) 
in the domain where < r < /3, < r' < (3. Keeping r' fixed, then 

G{x, 0; x', r') = -|rr(e-^^V+(^', r')iP{^, 0)) 

Zj 



i 
Z 

Q_, 

qG{x,P;x',r') . (4.12 a) 



= -|Tr(V'(f,0)e-^^V'+(^',r')) 
= -|rr(e-^^V^(f,/3)V^+(f',r')) 



Similarly, 

G{x, r; f , P) = qG{x, r; f , 0) . (4.12 b) 

Since the hamiltonian is time- independent, G depends only on the difference t — t' and 
eqs. (4.12) may be written as 

G{x, x'; r - r' < 0) = qG{x, f ; r - r' + /5) . (4.13) 

Note that when g = l(g = — l),G'is periodic (antiperiodic) with period (5 in the range 
< r, r' < f3, which corresponds to the usual bosonic (fermionic) case. 

Now the next step is to obtain an appropriate Fourier decomposition in the variable 
T for the Green function satisfying the boundary conditions (4.13). For this purpose, we 
introduce a new function 

g{x,x',T)^q^/''G{x,x',T) , (4.14) 
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which according to (4.13), is periodic in r variable with period /3. Thus this function can 
be expanded in a Fourier series 

g{x, x', r) = W e-'^-^g{x, x' , , (4.15) 

P n 

where cun — 2mr/p. The associated Fourier coefficient is then given by 

g{x, x', a;„) ^ J dr e'^-^g{x, x', r) . (4.16) 



Eqs. (4.15) and (4.16) can be written in terms of G, 

G{x, x', r) = W e-''^^g{x, x', h) , (4.17) 

P n 

g{x, f , h) = Jdr e'^^^G{x, x\ r) , (4.18) 



where 

A;4 = 2n7r//3- ^Ing . (4.19) 

P 

When q = I {q = —1), then A;4 = 2nn / (3 {k^ = (2n + l)7r//?), i.e., we have the common 
expressions for the bosonic (fermionic) frequencies. 

For the noninteracting Green function (4.9) we have then 



g°{x - x', h) = Jdr e'^^^G°{x - x\ r) 



1 - ? 



k 

= ^ E {{-'^'''~"'^' + l) < ^^^^ >} ■ (4.20) 

k ^ 

The condition (4.19) and eq. (4.11) imply that the expression in brackets (...) in (4.20) 
is equal to 1 and consequently, we have 

^ik{x—x') 

9%x-x',h) = -Y.j^-^^, 

k 

or after performing the Fourier decomposition in the space variables, 

g°(k,h) = —^. (4.21) 

We see that the free g-propagator has the usual form with the only difference that now 
the frequencies are given by eq. (4.19). 
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To conclude this section let us discuss the method for evaluating frequency sums 
appearing in the calculation of statistical averages. Suppose that for a function F(..., /C4) 
we want to evaluate the sum ^ ^F{..., k^). By choosing a function having poles at the 

values z — ki, and with unit residue we can reduce the summation over k^ to evaluation 
of a contour integral. One possible choice are the functions 

Then we have 

-^Y.F{...M) = -^E/<f(^.) , (4.23) 

where Zj are the poles of F in the complex plane defined by the variable z — k^^. Let 
us remark also that the choice of or is determined by the requirement that the 
integral / fFdz should converge on the circumference of an infinite-radius circle in the 
complex plane. 

As an example of the use of formula (4.23), let us calculate again the thermal average 
< A'^Ak > but this time by means of the Green function. 

Prom the definition (4.1) and the expansions (4.8), it is easy to show that 

< At A, -- - Y: , (4.24) 

where e > 0. (The factor e*'^"^ is a convergence factor and therefore e must be kept until 
the sum is evaluated). Using the formula (4.23) with the function defined in (4.22), 
we get 

1 



of&k - q 



i.e., we reproduce the same result obtained in the previous section (see eqs. (3.9) and 
(3.11)) from statistical mechanics. 

To develop an analogous Green function method for anyonic fields is of considerable 
interest. This question is under study. 



5 The ideal q-gas 

By an "ideal g-gas" we understand a system defined by the hamiltonian 

H ^^ei {a aftti + (1 - a) aiaf) , (5.1) 

i 

where the operators Oj, af satisfy the g-oscillator algebra (2.21), and a is a real parameter 
between and 1. An equivalent form of the hamiltonian can be written using (2.4): 

H^Y.^i + (1 - + 1]) ■ (5-2) 
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We shall interpret ai,af,Ni as annihilation, creation and occupation number operators, 
respectively, of particles in the state (level) i, although the mathematical results are, of 
course, independent of this interpretation. We shall cal the energy of the level i. The 
ideal q-gas is a special case of a general class of systems, which could be called "general 
ideal gases". The energy eigenvalues of such a system are uniquely determined by a set 
of occupation numbers {n,} and a set of functions Ef. 

E = Y. Hni) (5.3) 

i 

In a usual ideal gas, the functions Ei are proportional to the occupation number: 

Ei{n) = e°n , (5.4) 

corresponding to the physical interpretation of the total energy as a sum of single-particle 
energies e°, with vanishing interaction energy between the particles. In the general case 
(5.3) the interaction between the particles is such that it causes the energy of the ith level 
to depend on the occupancy of that level. For the ideal g-gas (5.2), this dependence is for 
real q = exp t 

EAn) = -^(a sinh (nt) + (1 - a)sinh ((n + l)t)), (5.5) 
sinht 



and for q = exp{ii 



Ei{n) = -^{asm{n9) + (1 - a) sin((n + 1)9)) . (5.6) 
sin 9 



(In an interpretation, where Oj, af and Ni are the ladder operators and excitation levels 
of oscillators, (5.5) and (5.6) correspond to a certain type of nonlinear oscillator |]25|). 
The grand canonical partition function 

Z = Trexp{-(3{H - fiN)) = exp(-/3fi) , (5.7) 

where is the total number operator 

N = J2N^, (5.8) 

i 

/i the corresponding chemical potential, and Q the grand canonical potential, factorizes 
for a general ideal gas into a product of single level partition functions: 

Z = llZ,{z,(3,fi) , (5.9) 

i 

CO 

Z^{i, p,fj.) = Y. e-^(^^(")-'^") . (5.10) 

n=0 

Correspondingly, the grand canonical potential is given as a sum over the levels i: 

^ = -^T.^ogZ,{z,(3,fi) . (5.11) 

The allowed values of the chemical potential fi are determined by the requirement that 
the sum in (5.10) converges. For the ideal g-gas with real q, the absolute value of Ei grows 
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exponentially with n, eq. (5.5). In order to have the total energy bounded from below, 
we require > for all i in (5.5). Then fj, can take any real value if all > 0, if at least 
one ej = 0, we must have /j, < 0. When q is a pure phase, |-E'i(7T')| stays bounded for all in 
(eq. (5.6)), in this case we require fi < 0. 

The sum over n in (5.10) cannot be explicitly performed in general. However, when 
q = expt is close to 1 (t — > 0), we can calculate the leading correction to the usual Bose 
gas. To order t^, eq. (5.5) reads 

Ei = ei{n + l-a + ^n{n + l)(n + 2 - 3a)t'^ H ) . (5.12) 

Substituting this into (5.10) gives the approximate expression 

Zi{i, (3, n) ~ e^("-^)^' e^('^-'')"(l - -I3ein{n + l)(n + 2 - 3a) + • • •) ■ (5.13) 

n=0 6 

Introducing the variable Xi = (3ei and the fugacity z = exp(/3/i,) and performing the sum 
over n we obtain 

Zi{xi,Z) ~ e(°-^)^'Zo(^,Xi)(l - t'^XiZe-'''Zl{z,Xi){aze-''' + 1 - a)) , (5.14) 
where ^ 

M^,^) = — 

1 — ze ^ 

corresponds to the ideal Bose gas. The grand canonical potential is then to second order 
in t 

~ 5^(1 - a)xi - ^ log Zo{z, Xi) (5.15) 

i i 

+t^zYxie~'''Zl{z,Xi){aze~''' + 1 - a) . 

i 

In order to study the thermodynamics of the system we have to specify the energies 
and the density of states. We shall consider a system where the states i are specified by 
a d-dimensional momentum vector k (the momentum of the g-boson) . The energy of the 
g'-boson is given by the dispersion law 

ei = e{k) = ^\kW (5.16) 

covering the cases of nonrelativistic (7 = l/2m, p = 2) and ultrarelativistic (7 = 1, p = 1) 
g-bosons. We enclose the system in a large d-dimensional volume V and replace in the 
usual way the sum over levels by an integral over /c-space: 



{2nY 



The first term on the right hand side of (5.15) is now divergent, giving rise to an infinite 
(negative) vacuum pressure, and we renormalize it away by simply dropping it. This 
means in effect that we subtract away the vacuum energy, i.e. replace the functions Ei{n) 
by Ei{n) - Ei{Q). 
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Performing the fc-space integrals using the general formula 

I 



oo a;"-ie-^'^ T(a) ^ (k + n - 1)\ z'' 

dx — r- = 7TT 2^ 



{el'' - zY (n-1)!^ k\ {qk + qn + pY 
and introducing the functions ge{z) 



(5.17) 



oo 



k=l 



(5.18) 



we get for the pressure p— —il/V 



Pp=aLi^,{z)--^t^ 



gd_i{z) + (1 - 2a)gd{z) 



Here 



_ Aa_, nj + 1) 
where A^-i is the area of the unit sphere 5''^"^: 




d/p 



(5.19) 



(5.20) 



27^ 

r(i) 



The g'-boson density 



is easily found, since 



n 



N 
V 




Ty 



(5.21) 



dgt{z) 
dz 



We get 



n — a 




g<L. 



+ (1 - 2a)gc 



(5.22) 



According to our general considerations, the chemical potential is restricted to negative 

values for the system we consider. 

Since 

ge{z) > ge'{z) 

when £ < i' and z ^ 1, eq. (5.19) implies that the pressure becomes negative for n 
sufficiently close to {z sufficiently close to 1) if t is real. This is not a signal of instabihty 
of the ideal g-gas, rather it signals the breakdown of the expansion in t. From (5.13) it is 
clear that the effective expansion parameter is times a positive power of ZQ{z,Xi), and 
this grows large when z ^ 1 and — > 0. Although the problem of negative pressure is 
absent for imaginary t, also in this case the expansion cannot be trusted as z — > 1. Thus 
we are not able to address the interesting question how the phenomenon of Bose-Einstein 
condensation is affected by the deformation. 

For low densities, however, the expansion in is meaningful. We invert eq. (5.22) to 
obtain z as a series in n: 

Tl 71 

z ^ a,{-) + a2{-)^ + ■ ■ ■ , (5.23) 
d a 
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where 

ai = 1 + -(1 - , 
P 

a = T^(l at ) . 

2d/p ^ p ' 

Substituting (5.23) into (5.19) and expanding in powers of n, we obtain the virial expan- 
sion of the equation of state 

(3p = n{l + Bn + ---) . (5.24) 
The second virial coefficient is given by 

B = — -i_(i--t2N (5 25) 

Thus, for very low densities, the g-gas behaves as a classical ideal gas. The expression for 
the second virial coefficient is independent of the parameter a in the hamiltonian, and 
shows that the deformation weakens the attraction between pairs of g-bosons producing 
an increase in the pressure, when t is real, whereas the effect of the deformation is the 
opposite, when t is imaginary. 

Our results in this section generalize previous partial results obtained e.g. in refs. 



|2J] and [^. The papers in ref. address the same questions as we have done in 
this section using the retarded and advanced Green functions. Their results, where both 
the chemical potential ^ and the particle number N enter into the average distribution 
number and have sharp values simultaneously are meaningless, however. 



6 Anyonic gas near Bose and Fermi statistics. 

Some of the statistical properties of the anyonic gas have been studied in certain ap- 
proximations In this section we will find perturbatively (in the first order of 
the statistical parameter) the equation of state for a gas of anyons near Bose and Fermi 
statistics by using instead the temperature Green function method. The same results 
have been obtained in [BO] in a different approach. 



We consider a system of N identical anyons described by the hamiltonian |15 



N 



{Pi - Aj) 
2m 



where 



A,, 



r 



(6.1) 



(6.2) 



is the statistical gauge field, k is a unit vector perpendicular to the plane and a = e^/vr 
is the statistical parameter (e and are the charge and the fiux carried by each anyon). 

Due to the singular a^-interaction term in (6.1), the standard perturbation treatment 
requires a redefinition of the A^-body anyonic wave function, so that the latter vanishes 



when any two anyons approach each other |29]. Therefore, the real wave function is 
defined as 



n '^ij'v'p(^i,--,nv) 



i<j 
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For the new function ipp the hamiltonian reads as 



N ( -^2 ■ U ~' II 



=1 [ 



2m m rfj m ^ rf^ 



(6.3) 



and consequently, the singular terms diverging like ^/rfj are cancelled. Now the pertur- 
bative analysis can be carried out. At order a, only the term proportional to |a| in (6.3) 
contributes to the grand partition function. Indeed, the contribution of the second term 
proportional to a can be shown to vanish due to the symmetry of the free spectrum under 
the change px Py (i.e., the corrections to the spectrum coming from the a interaction 
cancel in the order a in the grand partition function). Furthermore, the |a| interaction 
term in (6.3) can be replaced (in a perturbative sense) by a sum of two-body 5-potentials 



29 



y^2j[Mj26\f,-f,) . (6.4) 

^ i<j 

This fact will allow us to simplify the calculations of the self-energy corrections. 

Let us consider now the statistics. As is well known the temperature Green function 
satisfies the Dyson equation, which in the momentum space has the form |^ 

G{k,Un) = G\k,Un)+G\k,UJr,)Y.{k,Un)G{k,UJr,) , (6.5) 

where Yl, is the proper self-energy and G° is the free propagator, 

G° = ^ , (6.6) 

ojn = 2?27r//9 {ujn = (2?2+ l)7r//3) for bosons (fermions), G^= P/(2m) and fj, is the chemical 
potential. Eq. (6.5) has the explicit solution 

G{k,UJn) = ^ . (6.7) 

iuJn - (Sfc -A*) - Eik^UJn) 

The first-order self-energy is easily evaluated and is given by 

d^k 



±/7^n^-^'K', (6.8) 



{2n) 

where the upper (lower) sign corresponds to bosons (fermions), 

nl = N (6-9) 

" e^(^fc"^) =F 1 

is the usual Bose (Fermi) distribution and V is the two-particle interaction potential. 
According to (6.4), for anyons we have 

V(k) = —\a\ (6.10) 
m 
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and therefore, J2{i) is easily evaluated. Finally we obtain for bosons 



sf« = yy<'- — =-^'"(1-.). (6.11; 



2\a\ r , z 1\ol\ 



Sfi) = . (6.12) 



while for fermions, 

We conclude that there are no connections of order a to the thermodynamical quantities 
of the anyonic gas near Fermi statistics. 

Since Yl^\) is independent of the particle distribution will be given by the same 
expression as the unperturbed one (6.9) with the only difference that instead of the free 
energy G^, we will have +Z]^)- Besides, in the approximation we are dealing with 

(first order in a), Yl^\) is also independent of k and thus, the density of particles is easily 
found. 



n 



f d^k 1 , -aV^ ^ 



1 o 1 f 2/3 Em 

- -y, ln(l - z + z/5Sf,)) ~ -- |ln(l - z) + 

where A = . 

V m 

Substituting eq. (6.11) into the latter equation, we obtain 

n = -lln(l-.)|l-^| . (6.13) 



A 

For the pressure p, in turn, we find 



Pf3= I dz- = -\Y.^-\^\ -^) \ ■ (6.14) 



Eqs. (6.13) and (6.14) were found in p9| by following a different procedure. These 
equations can be rewritten in terms of the functions gi{z) defined in (5-18). We have then 

n=^g^{z){l-2\a\go{z)} , (6.15) 

pP = ^{92{z) - \a\gl{z)} . (6.16) 
Expanding p in powers of n we get the following expression for the second virial coefficient , 

B = + |a|A2. (6.17) 

This expression shows that B exhibits a nonanalytic behaviour (a cusp at a = 0) as a 
function of the statistics determining parameter a. 

Comparing eqs. (5.19), (5.22) and (5.25) for the two-dimensional non-relativistic g-gas 
with eqs. (6.15), (6.16) and (6.17) for the anyonic gas, we notice that they differ from 
each other and consequently, one sees explicitly that these two objects indeed describe 
different physical systems. 
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7 Conclusions. 



In this paper wc have studied some statistical properties of g-oscillators, quons and anyons 
in order to shed some hght on the interrelation as well as on the basic differences among 
these objects. In particular, the statistics of bosonic and fermionic g-oscillators have been 
considered in detail for different Fock space representations of the corresponding algebras 
in the case when the hamiltonian is identified with the number operator. The choice of 
this hamiltonian allowed us to calculate the statistical averages and also to introduce the 
Heisenberg picture of operators. The Green function method has been also introduced for 
nonrelativistic g-particles by assuming that Green function satisfies g-deformed periodicity 
conditions. This general method is hoped to be employed in attacking the problem of 
formulating the statistics of an anyonic gas. 
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